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We study the photonic interactions between two distant atoms which are coupled by an optical 
element (a lens or an optical fiber) focussing part of their emitted radiation onto each other. Two 
regimes are distinguished depending on the ratio between the radiative lifetime of the atomic excited 
state and the propagation time of a photon between the two atoms. In the two regimes, well below 
saturation the dynamics exhibit either typical features of a bad resonator, where the atoms act as 
the mirrors, or typical characteristics of dipole-dipole interaction. We study the coherence properties 
of the emitted light and show that it carries signatures of the multiple scattering processes between 
the atoms. The model predictions are compared with the experimental results in J. Eschner et at, 
Nature 413, 495 (2001). 



I. INTRODUCTION 

Control of photon-atom interaction lies at the heart of 
quantum technologies based on atomic and photonic sys- 
tems [H . Recent experiments demonstrated the quantum 
correlations between atoms and emitted photons 
Atom-photon entanglement was then applied for entan- 
gling distant atoms by photon measurement Q. Fur- 
ther experiments demonstrated the possibility to spa- 
tially confine atoms with nanometric precision inside res- 
onators [1, 0) [1] ) and hence to control their coupling with 
the electromagnetic field modes of cavities. Such preci- 
sion has permitted realizing quantum light sources with 
high degree of control d, [TO, [HI, [H, [S] , and hence to 
pose the basis for the realization of quantum networks 
based on atom-photon interfaces 

Parallel to these experimental efforts, studies are also 
focussing on achieving strong coupling between atoms 
and photons by means of optical elements, such as lenses 
of large numerical aperture [H, [l^ [H, \1T\. [isl . [l^ or 
optical fibers [l^, [2l|, |22| . In particular, in |14|| two dis- 
tant atoms in front of a mirror were coupled by means 
of a lens, focussing the radiation emitted by one atom 
onto the other. In this setup, the first-order coherence 
was experimentally studied, showing an interference pat- 
tern when the optical path length between the atoms was 
varied. In earlier experiments with two trapped ions, far- 
field interference of their scattered light [23'|, and their 
near-field interaction [23 | were studied. 

In this article we present an extensive theoretical study 
of the radiative properties of two distant atoms when 
they are coupled via an optical element, which could be 
an optical fiber or a lens, as sketched in Fig.[T] In this sit- 
uation radiation is multiply scattered between the atoms, 
until it is finally dissipated into the external modes of the 
electromagnetic field. Our model is based on the theory 
developed in [25,, J6] for the case of a single atom interact- 
ing with itself via a mirror, and extends it to the situation 
of two coupled atoms. The theoretical predictions of our 
model reproduce the experimental results of [14] and al- 
low us to identify possible measurements that highlight 
the multiple-scattering features. Moreover, the scattered 



photons are correlated with the scattering atoms, thereby 
establishing correlations and, in certain cases, entangle- 
ment between their internal excitations. 




FIG. 1: The dipolar transitions of two atoms, which are sev- 
eral optical wavelength apart, are coupled by a lens focussing 
part of their emitted radiation onto each other. In [3| a simi- 
lar situation was realized, coupling two atoms via a mirror and 
a lens. Analogous dynamics can be observed when the atoms 
are trapped close to an optical fiber, see for instance [2^ . 

This article is organized as follows. In Sec. |TT] we 
make some preliminary considerations on the system. In 
Sec. mil we introduce the model in detail and solve the 
basic equations describing the coupled dynamics of the 
internal atomic states and few photons of the electro- 
magnetic field. In Sec. IIVI we investigate in detail the 
radiative properties of the system, and in Sec. |V] we pro- 
vide the details of the first- and second-order coherence 
of the light scattered by the atoms when they are weakly 
driven by a laser. In Sec. I VII we provide some outlooks to 
the present work, and in the appendices we report details 
of the calculations. 



II. PRELIMINARY CONSIDERATIONS 

The scattering cross section of an atomic dipole transi- 
tion in fr ee sp ace is on the order of the square of its wave- 
length A |27| . Consequently, the free-space photonic in- 
teraction between two atomic dipoles at distance d is de- 
termined by the ratio X/d [1^: when d <C A, strong modi- 
fications of the atomic emission spectrum of one atom due 
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to the presence of another one are observable [2J, |29|, |30| ; 
when d X, these effects are neghgible, and the atoms 
scatter photons independently. This behavior is dramat- 
ically modified if an optical system, like a lens with large 
numerical aperture or an optical fiber, focusses a signifi- 
cant fraction of the radiation emitted by one atom onto 
the other. This latter situation is sketched in Fig. [1] for 
the case of a lens that images the atoms onto each other. 

When the photonic interaction between the atoms is 
mediated by an optical element, its strength is charac- 
terized by the fraction k of modes of the electromag- 
netic field which the optical system transforms into each 
other. Thus k replaces the scaling with X/d of the free- 
space case, and coupling over much larger distances than 
A may be achieved. 

The atom-atom distance d, or more precisely the prop- 
agation time for a photon from one atom to the other via 
the optical element. 



d 

T — — 

C 



(1) 



remains an important physical parameter of the photonic 
interaction, since it has to be compared with the radia- 
tive lifetime of the atomic dipolar transition I/7, which 
determines the time scale on which the photonic excita- 
tion is dissipated into free space, as well as the length 
of the emitted photonic wave packet. When 7r ^ 1, 
the process of photon scattering by each atom is well 
localized in time and space: a photonic excitation is ex- 
changed between the atoms at integer multiples of the 
delay time t, until its amplitude is damped to zero by 
emission into the external modes of the electromagnetic 
field. When jt ^ 1, in contrast, multiple scattering 
events add up coherently during the excitation time of 
each atom, causing the spontaneous emission rate to be 
enhanced or suppressed, depending on the interatomic 
distance (modulo the wavelength). This regime is equiv- 
alent to dipole-dipole interaction with a delay time r. 

In all cases, the system of two atoms confining radi- 
ation by multiple scattering shows some analogies with 
an optical resonator with low-reflectivity mirrors. This 
analogy is appropriate when the atomic transition is not 
saturated. Indeed, in this regime the radiative properties 
are very similar to those of a sing le atom interacting with 
itself via a mirror, studied in [l4 [Tsl. [26t . The pecuHarity 
of the two-atom system becomes more evident when sat- 
uration effects are relevant. Some important properties, 
such as the creation of correlations and entanglement be- 
tween the atoms via the multiply scattered photons, are 
identified when studying intensity-intensity correlation of 
the light scattered by the two-atom system, as discussed 
in Sec. IVl 



III. THE MODEL 

In this section we develop the theoretical model for 
describing the dynamics of two atoms in presence of an 



optical element which focusses the radiation emitted by 
each atom into the other, as sketched in Fig. In par- 
ticular, we use the theoretical formalism in [26| for one 
atom in front of a mirror, and generalize it to the case of 
two coupled atoms. 

The system consists of two identical atoms of mass M, 
which are trapped at the positions ri and r2 , and whose 
relevant electronic degrees of freedom are the ground 
state \g) and the excited state |e) forming a dipole transi- 
tion with dipole moment D, frequency uq and wavelength 
A — 27rc/wo- The interatomic distance d — |r2 ^ ri| is 
such that d 3> A, thus free-space dipole-dipole interaction 
between the atoms is negligible. We assume, however, 
that a lens (or an equivalent optical system) is placed 
between the atoms, which collects a fraction of the radi- 
ation from each atom and focusses it onto the other one. 
We use the plane wave decomposition for these modes 
and label them with p, in order to distinguish them from 
the external modes which do not couple the atoms; the 
latter are labeled with n, see also Fig. [5J The Hamilto- 
nian of the system describes the interaction between the 
dipoles and the modes of the electromagnetic field, and 
can be decomposed into the sum 



H = Hn 



Vr 



cmf : 



(2) 



where Hq gives the self-energy and Vomf the interaction 
between the dipoles and the modes of the electromagnetic 
field. In detail. 



Ho 



J = l,2 



E 



(3) 



where the first term describes the energy of the atoms, 
with (Tj = |.9)j(e| and crj its adjoint, and subscript 
j — 1,2 labeling the atom. The second term is the 
free Hamiltonian of the transverse photon field where the 
summation runs over all field modes. We label by £ the 
mode with wave vector and polarization ei -L k^, while 
aj and ai are the creation and annihilation operators for 
a photon in that mode, obeying the commutation rela- 



tion 



Of, al, 



Sg^i'. In particular, the modes with label 

£ = p are the ones which couple the atoms via the lens. 

The interaction of the atoms with the electromagnetic 
field, V^nif, is given in the electric dipole and rotating 
wave approximation, and takes the form 



E 



-l-H.c. 



(4) 



where 



= {D ■ ee) ^Lj,/{2eohV) 



with the vacuum electric permittivity Eq and the quanti- 
zation volume V. 

In presence of a laser driving the atoms the Hamilto- 
nian will be given by 



H'^H + VL{t)., 



(5) 
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where the term Vl describes the atom-laser couphng and 
reads 

VL^hnY, cr+e'(''^ ""J-"^*) + H.c. (6) 

Here, the laser is a classical field at frequency ujl |27i], 
is the coupling strength, and Icl is the wave vector of the 
incident laser beam. 

The dynamics of the system is studied by solving 
the Schrodinger equation treating the interaction of the 
atoms with the electromagnetic field as a perturbation. 
For this purpose, the wave function {ipit)) of the atoms 
and the field at time i, in the interaction picture with 
respect to Hq, is described by 

m))=b^'Ht)\ei, 92,0) +b^^\t)\g^, 62,0) (7) 
+ E 52, Ip, 0^) + b^^^ {t)\gi,g2. Op, 1^) , 

P M 

where the state |0) corresponds to the vacuum state of the 
electromagnetic field, and the state jn^^) {\np)) to n pho- 
tons in mode fj, (p). In Eq. ([7]) we have assumed that at 
most one excitation is present in the system. In particu- 
lar, the coefficients foi"'-' (t) are the probability amplitudes 
at time t for atom j being in the excited state, while the 
coefficient hf\t) gives the probability amplitude to find 
a photon in the field mode t at time t, with both atoms 
in the ground state. For later convenience, we also intro- 
duce the probability amplitudes bg''''^\t), with 

feW(<) = 6^^)(t) + 6f^)(t), 

and which distinguish which atom has emitted the pho- 
ton into mode £. 

We will solve the Schrodinger equation using this 
ansatz first in absence and then in presence of a laser 
driving the atoms. In particular, we will study the dy- 
namics as a function of two important physical quantities 
which characterize the system. The first is the time de- 
lay r for light to propagate from one atom to the other, 
defined in Eq. ([1]). As noted before, we consider the case 
CT 3> A. The second important quantity is the strength 
of the photonic coupling between the atoms mediated by 
the lens, which is defined through the fraction of Att solid 
angle within which the radiation from one atom is fo- 
cussed onto the other. This corresponds to the fraction 
of modes labeled with p, which propagate from one atom 
to the other via the lens. We denote the coupling by the 
dimensionless parameter k, 

n = E(l-|D.n,|V|Dp) 

^ A 1 df7o(l-|D.n|V|Dp) , (8) 

where Up = kp/fc and is the solid angle collected by 
the lens. The value of k lies in the interval < k < 1, 



whereby k — > corresponds to the limit without the lens 
and K — > 1 would describe an ideal optical system that 
maps all radiation from one atom onto the other. 

A. Perturbative solution of the Schrodinger 
equation in absence of the laser. 

When the atom-laser coupling is set to zero, then 
in the reference frame of the atoms the coefficients 
be ^ , b'^J'''^ , bg^'^^ obey the differential equations 

b(j){t) = -^gpe"'''"-Je'('^""'^'''*6^'')(*) 
p 

bi^'^Ht) = .9pe-*''"'Je-'('^°-'^'')*6(^')(i), (9b) 
b'g^''"Ht) = .9Me"*^'''Je-'("«-"^)*6(^')(i), (9c) 

where in the regime |r2 — ri| ^ A we have neglected 
processes in which a photon emitted into a mode fi by 
one atom is reabsorbed by the other one. 

A closed form for the coefficients of the dipole excita- 
tions is found by summing over the modes of the elec- 
tromagnetic field and by applying the Wigner-Weisskopf 
approximation as in [26|. The details of the calculation 
are reported in App.lXl The resulting equations take the 
form 

bi'Ht) = ~lbi'Ht)->^le'-"^b^^Ht-r)e{t-r), 

(10a) 

(10b) 

Equations pOap - (|10b[) show different behaviour depend- 
ing on whether t < t or t > t. For t < t these equations 
are decoupled and describe exponential damping at rate 
7 of the single-atom excited-state occupation, as in free 
space. After the time t, the coupling by light scatter- 
ing from each atom onto the other appears, its strength 
being set by the parameter k. 

We proceed by solving Eqs. (|10a[) - (jlObp for an arbi- 
trary initial state with a single atomic excitation, 

1^(0)) =ai|e,5,0)+a2|5,e,0). (11) 

A simple solution is then found by using the decom- 
position into symmetric and antisymmetric coefficients 

c±it), 

6W(i) = (C+(i) + C_(t))/2, (12) 
6(2)(i) = (C+(t)-C_(t))/2, (13) 

obeying the differential equations 
C±(t) = ~lc±{t) T «|e'"«^C±(t - T)e(t - r) , (14) 
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whose solution is [45j 

oo 

C±W = C±(0)^(±l)^-4(t), 

with 



fc=0 



I,{t) = ^ ^ {t - krfe-i^'-^-)Q{t kr) . (15) 

Correspondingly, the probability amplitudes for the ex- 
cited states are 

b^^\t) = ai^hkit) + a2^hk+i{t) , (16a) 

k k 

bf\t) = ai^/2fc+i(i)+a2^/2fc(i), (16b) 

k k 

while the probability amplitudes bg-''^\t) for the emission 
of a photon into mode fi by atom j are given by 

,-ik„ri 



(17) 



^ [aiH2k{t,uj^) + a2H2k+i{t,u}f,)] 



k=0 



i+i^. 



(18) 



[aii/2fe+i(i, ^m) + "2^^2fe(t,t^M)] 



fe=0 



with 



(19) 



In Eqs. ([T7)) - (|18p we assumed that the electromagnetic 
field is initially in the vacuum state, 6^(0) — 0, and we 
introduced the function 



Hkit,Lu) 



(-Kfe'"^)'= 
A?! 



with tk = t — kr and 

Gk{s) = iFi(fc,fc + l,-s)-e- 



(20) 



(21) 



where iFijk, k + 1,— s) is the confluent hypergeometric 
function [47|. In the limit k — > 0, i.e. when there is 
no couphng between the atoms, Eqs. pT|) - P^ reduce to 
the usual free space deca y sp ectrum of two independent 
dipoles with linewidth 7 [46l |. 



B. Perturbative solution of the Schrodinger 
equation in presence of the laser. 

We consider now the situation that the atoms are 
weakly driven by a laser at intensity fi. Hence, we set 
n ^ in the Schrodinger equation and solve the dynam- 
ics of the new Hamiltonian assuming that Vl is a weak 
perturbation to the atomic dynamics. We use the ansatz 
for the wave function in Eq. ([7]), where we denote now 

the probability amplitudes by bi^\t),b\j''^\t),b^g\t) 
ci^\t),c\j^-'^\t),c^g\t) (with j = 1,2 and t = n,p). Let 
IV'(O)) = 1517 52,0) be the initial state. By solving the 
coupled differential equations for the probability ampli- 
tudes in first order in and in the reference frame rotat- 
ing at the laser frequency ojl we find 



«(t) = -i/a / d<'e--*(ei,ff2,0|e-'^(*-*')/''Fi(i')e-'^*'/''|5i,52,0) 



/ dtV--(*-*''(ei,52,0| [e-'^(*-*')A (e*i--i |ei, 52, 0) + e*-'-^ 62, 0)) /\/2" 



(22) 



Corresponding expressions are derived for Ce {t) and can write 
c"p (t) . The term inside the square bracket corresponds to 
the time evolution of the state |/3o) = (e'*'^ '"^ |ei, ^2, 0) + 
gikL r2 g^, 0))/\/2 when there is no laser. Hence we 



c. 





t 

(^)(0 = -iV2n f dt'e-'^'''b[^\t') , (23b) 
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where f. = p, ^ and we have introduced the detunings 



(,0') or,^ ?,(^) 



(24a) 
(24b) 



The coefficients be and bg ' are found using the solutions 
derived in Sec. IIII AI when the initial state is \/3o)- One 
gets 



-ifl 



(25) 



J2 [H2kit,LUL) + e"^^ H2k+iit,u;L)] 



k=0 



with 



(/fL = kL • (r2 - ri) . 



(26) 



The equation for Ce (t) results from Eq. ((25| by inter- 
changing the indices 1 <-> 2. The probability amplitudes 
c'g'^\t) are found using Eq. p?)) in Eq. (|23b|) . assuming 
that initially both atoms are in the ground state and the 
electromagnetic field in the vacuum state. One gets 



di'e'("-c-'^f<)*'e'(''L-kM)-ri 



not decay appreciably. This shows in more detail how 
the radiative properties of the system are determined by 
the parameter 7r, the ratio between the delay time and 
the excited state lifetime. In particular, for 7T ^ 1 in- 
terference plays no role, and the photonic excitation is a 
wave packet bouncing between the two atoms, until its 
intensity is damped to zero by scattering into free space. 
For 7r ^ 1 the terms in psp add up coherently and in- 
terfere. The effect of the interaction hence modifies the 
radiative properties of the atoms, and the dynamics are 
analogous to an effective dipole-dipole interaction [2^. 

In this perspective, the optical set-up composed by two 
atoms and the lens can be considered like a resonator, 
where the atoms arc mirrors of low reflectivity and re- 
flection bandwidth 7/2, while 2t is the round-trip time. 
The parameter 7r hence gives the number of modes that 
this peculiar "two-atom cavity" sustains: for 7r 3> 1 it 
sustains several modes and can be considered a "multi- 
mode resonator" . Conversely, for 7r <C 1 only a single 
mode of radiation is supported, and we will denote this 
case as " single-mode resonator" . 

Using this insight, we now analyze the probability am- 
plitude and the spectrum of the emitted photons in the 
external modes labeled with /i. Let us first assume that 
the laser is absent, and that initially atom 1 is in the 
excited state, i.e. ai = 1,Q!2 = in Eq. pT|) . From 
Eqs. (fT7)) - p8|) . the amplitude probabihty for the state of 
the field reads 



k=0 



with (5^ given in Eq. (jl9p . The probability amplitude 
cf'^\t) for atom 2 is obtained by swapping the super- 



scripts 1 ^ 2 in Eg. (|2^ 



Discussion 



The probability amplitudes of the atomic excited states 
in absence and in presence of the laser, given in Eqs. (I16|) 
and (|25p , respectively, are the coherent sums over contri- 
butions starting at different instants of time — t ~ kr. 
These contributions correspond to the effect of k ex- 
changes of a photonic excitation between the two atoms. 
In particular, for the case of atom 1, the contributions 
at T2k correspond to an excitation which propagated to 
atom 2 and back. Hence, in Eq. (fT6|) this term vanishes 
when initially only atom 2 is excited. Similarly, the con- 
tributions ai t = T2k+i vanish when atom 2 is initially 
in the ground state. Similar considerations apply for the 
case in which the laser drives the atom, Eq. (|25p . 

An important property of these equations is that each 
term of the sum has a well-defined phase, which is an 
integer multiple of iOgT {ujlt with the laser excitation). 
At the same time the contributions are damped by an 
exponential function at rate 7. Consequently, the indi- 
vidual terms show interference if over the time r they do 



5m 



■E 



(28) 



fc=0 



X e 



'H2k{t,UJ^,)+C 



-iku 



2k+l 



where the two terms under the sum account for the re- 
spective contributions of the two atoms to the emission 
into the field mode. The label k gives the number of pho- 
ton exchanges between the two atoms before the photon 
is finally emitted into the external mode /i. 

In the long time limit Eq. (pS)) reduces to the form 



b^^Ht 



(1 +\S ) - ^-7piw^,r(l+cosi?) 

00)- a c-"^"-'-i ^2 ^ 

' ^'^ (i+i<5^)2-(K§e---)2 



where 'd denotes the angle between the vector ri — r2 and 
the wave vector of the mode, see Fig. [21 At = 7r/2, in 
particular, the probability to measure a photon in mode 
fi is given by 

|5('^)(i-oo)|2 = 



^[1 -f KC0SWpT]2 + [too 



sinaj„rl 



(29) 

showing that the spectrum exhibits a modulation at mul- 
tiples of the frequency 1/2t. In the resonator picture, 
the modulation peaks are at the mode frequencies of 
the resonator, and 1/2t corresponds to the free spec- 
tral range. The spectral modulation will be visible when 
7T ^ 1, i.e. when the system is in the "multi-mode- 
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resonator" regime. On the other hand, in the "single- 
mode" regime jt <^ 1, one will observe a change of the 
radiative linewidth, which depends on the phase luqt. 

When the atoms are laser-driven, the probability am- 
plitude for the excited state occupation in the long-time 
limit is 



.(1) 



(t oo) 



-in 



1 - Xe"^^ 



(30) 



while the probability amplitude that mode fi is occupied 
by one photon scattered by atom 1 takes the form 



X e 



i(kL 



5m 



2 

(1-Xe''^^) -fO(l) . (31) 



Here 



K = K- 



2i+iA 



(32) 



1 sin(At/2) 
A 



is the diffraction function [21 



and (5(*)(A) 

The second term on the right-hand side gives no contri- 
bution to the rate of emission, so that it is not explicitly 
reported. Its specific form can be found from Eq. ([77]) 
and Appendix |B] The probability amplitudes for the sec- 
ond atom are obtained by swapping the indices 1 <-> 2. 
The detailed derivation of these expressions is reported 
in Appendix [Bl 

These results are discussed for various specific limits 
in the following sections. 



IV. RADIATIVE PROPERTIES 

In this section we discuss the radiative properties of 
the two atoms, when they are observed together or indi- 
vidually, in the absence of laser excitation, and assuming 
that atom 1 is initially excited. The various quantities 
which will be discussed correspond to measurements with 
different detectors, as illustrated in the detailed set-up in 
Fig.H 



A. "Multi- mode- resonator" regime 

When 7r ^ 1, then the transient dynamics of the 
system is characterized by the two atoms exchanging a 
photonic excitation well localized in time. The excita- 



tion probabilities Pj = \bi''\t)\'^, with be^'{t) given by 
Eqs. p6)) . are displayed in Fig. [3] as a function of time. 
One clearly sees that a photonic excitation propagates 
back and forth between the atoms, while its amplitude 
is damped due to the scattering into the external modes 
of the electromagnetic field. The shape of the photonic 
wave packet exchanged between the two atoms changes 
with time: with each bounce it acquires a more symmet- 
ric and broader shape, due to the frequency-dependent 



Atom 2 




FIG. 2: Detailed schematic of the physical system showing 
the detectors which correspond to the various measurements 
described in the text. Another detector 2, /i would be placed 
in a location equivalent to that of detector to measure 
the emission of atom 2 individually. 



reflection by the atoms. The broadened wave packets 
increasingly overlap with time, such that interference be- 
tween subsequent excitations may become visible for long 
times, as shown in the example of Fig. 31 

The effects of the coherent addition of the multiple 
scattering events become more visible by inspecting the 
time-dependent probability of emitting the photon into 
the external modes of the electromagnetic field. In the 
continuum limit of Eq. (pS)) , it takes the form ^] 



(33) 



which for t — cxd coincides with the emission spectrum. 
An example of S{uj,t) is shown in Fig. O For times 
t < T, before scattering events can interfere, it exhibits 
a Lorentzian form like an atom in free space, while after 
a time i > r it develops spectral modulation with peaks 
spaced by l/2r. 

The effect of the distance between the atoms on their 
individual emission spectra is displayed in Figs. [HKa) 
and (b). In particular, the maxima of the spectra, spaced 
by the "free spectral range" 1/2t, shift according the op- 
tical distance between the atoms. The visibility of mod- 
ulation is larger, the closer k is to unity. 



B. "Single- mode- resonator" regime 

We now analyze the regime 7r <C 1, in which several 
photon excitations are exchanged between the atoms dur- 
ing the natural lifetime of the excited state. In Fig. [7]the 
excited state populations of both atoms are displayed. 
As atom 1 is initially excited, atom 2 stays in the ground 
state until the instant t = r, after which its excited state 
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FIG. 3: Excited state occupation of the atoms as a function 
of time in units of r, as given in Eq. (|16p . when initially atom 
1 is excited. The other parameters are 7r = 10, «: = 0.4, 
and ujoT = nvr. Note the change of vertical scale from each 
maximum to the next one. 




t/x 



FIG. 4: (Color online) Logarithmic plot of the excited state 
occupation of atom 1 as a function of time in units of r, 
as given from Eq. (|16p . The curves correspond to the values 
UJOT = nn (red solid line) and ujot = (2n+l/2)7r (green dotted 
line) for k = 0.4 and jt = 5. The dashed line corresponds to 
K = and is plotted for reference. 
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FIG. 5: (Color online) Probability of emission of a photon, 
Eq. (|33|l . as a function of frequency (in units of 7) and of time 
(in units of r), for the emission angle i? = -I and the phase 
uioT = 2n7r. 
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occupation increases due to the interaction with the radi- 
ation from atom 1, see FiglTfb). The excitation of atom 1 
is damped like in free space until time t = 2r, after which 
the damping rate is attenuated or enhanced depending on 
the relative interatomic distance, i.e. on ojqT. The effect 
of the relative phase between the multiple absorption- 
emission events is more evident when plotting the exci- 



FIG. 6: (Color online) Spectrum of the light emitted by atom 
j, Sjlco) cx limt-,ao \bg''\uj,t)f , as a function of the frequency 
(in units of tt/t) for (a) ujqt = nn and (b) uot — (2n + 1)-|. 
The parameters are k = 0.4 and 7r = 10. The dashed blue 
line gives the spectrum of the atom when k. — and is plotted 
for comparison. 



tation probabilities on a logarithmic scale, as shown in 
Fig.li 
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FIG. 7: (Color online) Excited state occupation of (a) atom 1 
and (b) atom 2 as a function of time (in units of r), as given in 
Eq. (|16p . when initially atom 1 is excited, and for ujot — nir 
(red solid line), ujqt = {2n + 1)-| (green dotted line). The 
dashed blue line is the solution for n — 0) and is displayed for 
reference. The other parameters are yr = 0.4, k — 0.4. The 
red-solid and green-dotted lines in the inset of (a) display the 
difference between the value of |&e^-'(t)p at ujqt = mv and at 
LJoT = (2?i+l)-|, respectively, from the corresponding excited 
state occupation at k = as a function of time. 



Figure [Hl^a) displays the emission probability S{u!^t) 
as a function of frequency and time, showing that it is 
always a single-peaked curve, whose width varies with 
time. Figure [Hb) displays the emission spectrum of the 
first atom in comparison with the one in free space for 
different values of the parameter wqt, showing that de- 
pending on the relative distance one can observe subradi- 
ant or superradiant emission. The atomic interaction is 
hence a retarded dipole-dipole interaction, mediated by 
the photonic excitation over the interatomic distance. 



C. Two atoms vs. single atom 

The cases studied so far share several analogies with 
the radiative properties of a single a tom in front of a mir- 
ror, analyzed for instance in (25l.l26||. In particular, in psl ] 
Alber studied the dynamics of one photon coupled to one 
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FIG. 8: (Color online) Logarithmic plot of the excited state 
occupation of atom 1 (solid line) and atom 2 (dotted line) 
as a function of time (in units of r) for the initial state 
1*^(0)) = 9, 0). The other parameters are 7r = 0.4 , n = 0.4 
and LuoT = (2n -|- l/2)7r. The blue dashed line shows the cor- 
responding atomic excitation for k = and is plotted for 
reference. 
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FIG. 9: (Color online) (a) Probability of emission of a photon 
by atom 1, Si{uj,t) oc \bg^\uj , t)\'^ , as a function of the fre- 
quency of the emitted photon (in units of 7) and of time (in 
units of t) for uiqt — nvr. (b) Spectrum of emission for atom 
1, Si{ul>), for ujqt = TiTT (red solid line) and loqt — (2n -I- 1)^ 
(green solid line). The other parameters are 7r = 0.4 and 
K — 0.4. The dashed blue line gives the emission spectrum in 
free space and is plotted for comparison. 
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atom at the center of a spherically symmetric cavity with 
perfect reflectivity. Depending on the radius of the cavity 
mirror, and thus on the time the photon needs to travel to 
the mirror and back in relation to the atomic decay time, 
Albcr defines the small- and large-cavity limit, whereby 
in the first case the atom-cavity system is characterized 
by a delocalized excitation, while in the second case a 
photonic wave packet propagates back and forth exciting 
periodically the atom. Although our system is a low- 
quality resonator, the multi-mode cavity that the two 
atoms form for 7T ^ 1 is analogous to the large-cavity 
limit in (H]. 

A very close analo gy e xists between our system and the 
system discussed in [26|, where Dorner and ZoUer inves- 
tigated the case of an atom interacting with its own light 
back-reflected by a distant mirror (26] . In particular, the 
dynamics of two atoms exchanging photons via the lens 
share strong analogies with the one of an atom inter- 
acting with its mirror image, if one restricts the Hilbert 
space to only one excitation, and if the atoms arc initially 
prepared in a symmetric state with ai = a2 = l/-\/2 in 
Eq. pT|) . For this initial state, the time evolution of 
the excitation of one of the atoms is the same as the 
one of the atom in front of the mirror. The probability 
amplitude for photon emission, however, shows some dif- 
ferences between the two cases. Figure [10] displays the 
emission spectrum as a function of the emission angle and 
of the frequency. Here, the oscillation of the intensity as 
a function of the angle of emission ?? is indeed a exclusive 
property of the two-atom case, arising from the fact that 
the light emitted from the two scatterers interferes in the 
far field. 



direction) and lolt (the optical path length between the 




(oOq - co)x/7t: 



V. LIGHT SCATTERING 

In this section we analyze the scattering properties of 
the system when the atoms are driven by a laser below 
saturation. In this case the time evolution of the ex- 
cited state amplitudes, Eqs. (|25p. describes the photon 
exchange between the two atoms, which now addition- 
ally interferes with the incident laser light. For 7r 3> 1, 
step-wise dynamics with the characteristic time step r 
are visible in the excited state occupation of each atom, 
as displayed in Fig. [TTT a) and (b). For different dis- 
tances between the atoms, and hence different phases 
of the various contributions, the discontinuities in the 
curves at multiples of r show constructive or destructive 
interference, while for long times t ^ t the excited-state 
population tends to a steady state value. In the limit 
7T < 1, displayed in Fig. [TlTc) and (d), the curves are 
smooth and tend to the same steady state values. This 
stationary value depends on the two phases tpL (the laser 



FIG. 10: Contour plots of the emission spectrum from both 
atoms S{uj), Eq. (|33|) . as a function of the frequency ui (in 
units of 7r/r) and of the angle of emission i9, for the initial 
state l'I'(O)) = -^{\e,g,0) + \g,e,0)). The other parameters 
are 7r = 10, «: = 0.4, and lvot = 2n7r (top panel), ujqt = 
(2n + 1)tt (bottom panel). 



atoms), according to 



|c«(i~.cx3)| = 



1-Ke 



|l_if2|2 ' ' • 

(34) 

where K is given in Eq. (|32p and is proportional to the 
coupling strength k between the two ions. Eq. ([M)) does 
not depend on the parameter 7r, which affects only the 
transient dynamics. When the atoms are not coupled, 
K = 0, one recovers the free-space steady state value, as 
found for an atom which is driven by a weak laser [27|. 



We note that for the specific value (pL — (2n -I- l)7r we which is the free-space formula with modified decay rate 
obtain and detuning, 

|cW(oo)|^^ , - , (35) 7 = 7(l-^cos.;ir), 

2 
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FIG. 11: (Color online) Excited state population of atom 1 when both atoms are driven by the laser, as evaluated from Eq. (I25|) . 
as a function time (in units of r). The figures are evaluated for k = 0.4, Q = 0.057, A = 0, and for = 20 (upper row) and 
for jT = 1 (lower row). The subplots (a) and (c) refer to the case ifL ~ {2n + l)7r, (b) and (d) to the case tfL = (2n — l/2)7r. 
The value of the phase ujlt for each curve is explicitly given in the plots. 



This result coincides with the excited state population 
of a single atom subject to interference between the laser 
excitation and the light back-scattered from a mirror [2^ . 
This equivalence holds only for the particular value = 
(2n + l)7r but not in the general case. 

In order to get some more insight, we analyze Eq. ([M]) 
for K <C 1. At first order in k the excited state population 
of the first atom takes the value 



.«(t-oo)| = 



72/4 + A2 
X (1 + 2kAcos [wlt + — 0]) 



(36) 



with 



A 



tan 




(37) 



The result for atom 2 is found from Eq. (|36p by swapping 
the subscripts 1 <-> 2, i.e. by changing the sign of ipL- 
Equation ([M)) shows how the excited state population is 
enhanced or suppressed as the parameter w^r is changed. 
This change in the atomic spontaneous emission rate, as 



well as a shift of the atomic resonance frequency, both 
controlled by the parameter lult, are manifestations of 
the modification of the radiative properties of the atoms 
due to their mutual interaction. Analogous frequency 
shifts in a single atom interacting with itself via a mirror 
have been experimentally observed by Wilson et al. p^ . 



A. Intensity of the scattered light 

Let us now consider the intensity of the light scattered 
by the laser-driven atoms, for several of the measurement 
set-ups illustrated in Fig. O 

First we consider the situation that the detection ap- 
paratus resolves the atomic positions and therefore sums 
up incoherently the photons emitted by the atoms (de- 
tector 1, /i plus detector 2, /i). Then the detection rate is 

= Ti,^ + Fa,^, whereby F,, ^ = limt-,^ \ci''''\t)\yt. 
Using Eq. ^ we find 



Ti.p = F2,p = 2TrS{uj^ - ujl) 



1 - Ke'^^ 



(21 

V 4 



|l-if2|2 



(38) 



where K is given in Eq. ([5^ . For K = 0, i.e. in absence 
of the optical element coupling the two atoms, the signal 
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reproduces the free-space resonance curve of the atomic 
dipole. For K it shows two modulations, with the 
phase 2llilt (through in the denominator) and with 
the phase lult + ip^ (in the numerator). The first one 
corresponds to previously emitted light returning to the 
same atom after scattering from the other one, the other 
modulation is produced by scattered laser light arriving 
from the other atom. In general, the modulations show 
how the scattering of a single atom is modified by the 
presence of another identical scatterer at a fixed distance. 
The maximum enhancement, when all scattering terms 
add up coherently, is found for 2ujlt — 27r, (p^ = and 
A = 0, and is equal to (1 + k)/(1 — k). For k = 0.2 it 
gives an enhancement of the signal of the order to 150%, 
as displayed in Fig. [121 For the case of a single atom inter- 
acting with itself via a distant mirror, analogous signals 
have been experimentally observed in Refs. p^. [Tsj. 



FIG. 12: (Color online) Intensity of the light emitted by one 
atom into free space, as a function of 2u!lt and for k = 0.2, 
when the system is driven by a laser at frequency lol = i^o and 
at ipL ~ 0. The intensity is normalized to the value obtained 
without coupling, k = (blue dashed line). 



^,,(K)/i;,,(0) 




When the light emitted by the two atoms is super- 
posed coherently on a detector (labeled fi' in Fig. [2]), the 
system is analogous to a double-slit set-up |23] with the 
important difference that the atoms additionally interact 
by photon exchange via the lens. In this case, the corre- 
sponding detection rate is T'^ = limt^oo \cg^\t)\'^/t and 
takes the explicit form 



r;, = SnSiuj^ - lol)-, — |cos[(kL - k^) • (ri - r2)/2] - Xcos[(kL + k^) • (n - v^)/2f . (39) 

I 



For K ~ 0, the observed spatial interference is the 
one of a double-slit set-up with two coherently driven 
sources [23., J33J. When K ^ 0, these properties are mod- 
ified by the multiple scattering. An important special 
case is when the direction of emission is k^ — — kL, cor- 
responding to the coherent backscattering direction , 
where one always finds a spatial maximum of the scat- 
tered intensity. 

We now consider a set-up in which one observes the 



modes p, through which the atoms interact (detectors 1, p 
and 2, p in Fig. [2]). This corresponds to the measurement 
arrangement in |l4!|. In this case, the rate at detector 

l,p is given by Fi^p = lim^^oo \c^g''Ht) + c[y \t)\'lt, 
superposing the light emitted by atom 1 directly into the 
detector with the light emitted by atom 2 towards atom 
1, and then into the detector, whereby kp and k^ are 
transformed into each other by the lens. It reads 



V^^p = Sn5{L0p-ujL)—, — |cos[(^L-WLr)/2]-ifcos[(^i+^LT)/2]|2 , (40) 



where we used that kp/ • r2 ~ kp • ri uj^t via the optical element. The corresponding rate r2.p is found 
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by changing the sign of tpL in Eq. (|40p . The total rate cal set-up, which was characterized by small values of n. 



AO)) 



.(0) 



^2°! '^as measured in Ref. [14] in an opti- Taking k ^ 1, the total rate T]^' reads 



(0) 



- 2,p 



r^*^-* cx 1 + cos cos wlt — kA [cos + 2 cos (/^l cos(cjlt — </))+ cos(2wlt — '/')] + O(k^) . 



(41) 



r 



where we omitted global constant factors, and and A 
are defined in Eq. p7|) . We observe that at zero order in 
K an interference pattern appears as a function of wlt, 
i.e. by changing the optical path between the ions. This is 
the classical interference of the light elastically scattered 
from both atoms into the same detector. The interference 
has visibility 



Vi = I cos ifL I 



(42) 



which is maximum when (pL = nn, with n integer, and 
which vanishes when (p^ = (2n + l)7r/2. This is a conse- 
quence of summing the signals from the two detectors, 
whose individual interference patterns may be shifted 
depending on ipL. It provides an explanation for the 
low contrast interference observed in the experiment of 
Ref. [l|. 

The vanishing contrast when the two signals are per- 
fectly anti-correlated provides a condition where the 
higher order effects in k, and thereby the interaction of 
the atoms, are particularly evident. Choosing this spe- 
cific condition, by varying the optical path length be- 
tween the ions one observes an interference pattern at 
twice the frequency of the classical interference, i.e. os- 
cillating with 2ujlt, with a shift determined by the 
detuning A, and whose visibility is given by 



V2 = kA 



^1 + 4A2/72 



(43) 



where we used Eq. (j37p . The visibility is maximum at 
atomic resonance. The doubled frequency of the inter- 
ference (compared to the classical one) with the inter- 
atomic distance shows that it is caused by two partial 
waves originating from the same atom, one reaching di- 
rectly the detector and the other being back-scattered 
once by the other atom. Analogously, processes where 
the same wave is scattered n times by the atoms give rise 
to interference terms with frequency uuilt and at higher 
order in k. 



B. Intensity-Intensity Correlations 



We now study the intensity-intensity correlations in 
this set-up, assuming that two detectors are placed in 
the far field of the scattered light at positions xi and 
X2 (corresponding to two detectors /i' in Fig. [2] at an- 
gles -di and -d^)- We denote by ^(2) (^i, t; X2, t + t') the 
(un-normalized) intensity-intensity correlation function 
for measuring a photon at time t and position xi, and 
another at X2 after an interval t' . It reads ,33] 



g(2)(xi,t;x2,t + t')- E e"=»"^^-"-")-*^+("^«-^<)-*=) (a+(i)4(t + i')'^.(^ + ^')'^cW). (44) 



where k — |kL| and x = x/|x|. Assuming that the atoms the two detection events. In this limit, we evaluate its 
are driven by the laser and have reached the steady state, explicit form in perturbation theory for the atom-photon 
Eq. (j44p depends solely on the time t' elapsed between interaction, and find the expression 
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(72+4A2)2 |1-X2|2 



(45) 



+ e' 



i(vi,-At') 



(1 — KcOSifL) 



where we have set 



ipj = fc(r2 — ri) • Xj = kdcosdj 



(46) 



and defined G^^^ {ipi, ip2;t') = g(2)(xi, xa, i + i')- The 
detailed derivation of Eq. (|45)) is reported in AppendixICl 
For K = 0, i.e. in absence of couphng, G'2) exhibits an in- 
terference pattern as a function of the distance |x2 ~ xi| 
between the detectors. Such interference emerges from 
two indistinguishable paths of two-photon emission. It 
has first been predicted in [sij for the case of two in- 
dependent quantum sources, and generalized in [s^l for 
the light scattered by two trapped atoms illuminated by 
a laser. In particular, the result of [s^l for weak laser 
intensity is obtained from Eq. ([45]) by taking the limit 
K ^ and setting A = 

We now consider this spatial interference pattern at 
t' = 0, for (fiL = and A = 0, but keeping k ^ 0. For 
these parameters it takes the form 



G(2)((pi,^2;0) = 



64^2^ cos2((^i - </'2)/2 



r 



1 



2k cosojo''" 



(47) 



In particular, Eq. ([i7|) vanishes for |v5i — iy92| = (2n-|- l)7r, 
showing strong antibunching at these points. Simi- 
larly, bunching is encountered whenever the condition 
It/Ji — '/32| = 2n7r is fulfilled. We note, moreover, that 
since this signal depends only on the difference ipi — if2, 
there exists a finite probability of measuring two pho- 
tons simultaneously at the positions of the screen which 
are the dark fringes of the first-order interference pat- 
tern, = {2n + l)7r. This behaviour has been discussed 
in [301 : it is connected to the fact that saturation effects 
diminish the contrast of the first-order correlation func- 
tion, leading to a non- vanishing probability of measuring 
a photon at these detector positions. The probability to 
measure the first photon in the dark fringe is essentially 
proportional to the occupation of the collective state 
|ei,e2), and the first detection projects the atoms into 
the antisymmetric Dicke state |V'(0)) = -^{\e,g) — \g,e)), 
which is an entangled state of the two distant atoms. 

The denominator of Eq. (|T7)l shows how the spatial 
interference pattern is modified due to the atom-atom 
interaction by multiple photon scattering, and how this 
modification depends on the phase lult. 

Figures [T51 and [Til displav the intensity-intensity corre- 
lation function versus t' and 1^92, for the situation 7r ^ 1. 
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FIG. 13: (Color online) G^^^ipi, ^52; t) as a function of time 
t (in units of r) and tp2 for tf \ = 2mi. The subplots refer to 
the cases (a) k = (no interaction) and (b) k = 0.4. The 
other parameters are Q, — O.O57, A = 0, = 0, 7r = 20, 
and uJoT = (2n -I- l)7r. 



The two figures correspond to the bright (tpi = 2mr) and 
dark {(pi = (2n+l)7r) fringes of the first-order correlation 
function, and both show the cases k = and k = 0.4, for 
comparison. One clearly observes the effect due to mul- 
tiply scattered photons, giving rise to abrupt changes 
in the slope of the correlation function at multiples of 
T. Hence, the interference due to multiple scattering en- 
hances or suppresses the probability of measuring the 
second photon at a certain time interval. Related effects 
have been observed in a single-atom interference exper- 
iment in Ref. (35| . We now analyze this latter property 
setting (fii = (2n -I- l)7r, i.e. when the first detector is set 
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FIG. 14: (Color online) Same as Fig.[T3]but with (fii = (2n + 
1)tt. The curves in (b) at t > r are magnified by a factor 30. 



at a dark fringe of the first-order correlation function: in 
Fig. [TlTa) one sees that for k = the second order cor- 
relation function is different from zero at = and for 
ip2 7^ 2n7r, and it vanishes after a transient time of the or- 
der 1/27, corresponding to the lifetime of the collective 
state |e,e) [S^. For k = 0.4, Fig. [MTb). one observes 
" revivals" when the time between the two detections is a 
multiple of T, and whose amplitude is strongly damped as 
a function of time. Inspecting Eq. (|45p for these specific 
parameters, we find that at short times the correlation 
function behaves as 



ipi = (2n+l)7r,t'<T 



(48) 



32rj4(i _ cos 932) 
7'*(1 + + 2kcosujot) 



and it is essentially proportional to the probability of 
measuring the atoms in the state \^{t')), obtained by 
freely evolving the initial state \ip{0)) = ■^{\e,g) — \g,e)), 
according to Eq. p^ . 

For longer times, t' > r, the second-order correlation 
function scales with and takes the form 



ipi = {2n+l)Tr,t'>T 

2 , 



(49) 



[1 + + 2kcosujot)\ ^ ^' ' 



which is essentially proportional to the stationary 
excited-state occupation of the atoms given in Eq. (PS)) . 

While the limit 7r ^ 1 is characterized by "revivals" 
of the correlation function versus the time t' between 
two photon detections, in the limit 7r < 1 one ob- 
serves a smooth decay of the correlation function with t', 
whereby its decay rate is modified depending on whether 
the multiply scattered waves interfere constructively or 
destructively at the atom. A comparison between the 
two regimes is shown in Fig. 1151 
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FIG. 15: (Color online) G^^^ ((^1, 1^2; t) as a function of time 
for i-Py = <4ii = (2n -|- l)7r and k = 0.4 when (a) 7r = 20 (note 
the change of vertical scale from each maximum to the next 
one) and (b) 7T = 0.4. The red-solid (bottom) and green- 
solid (top) curves are evaluated at ojot = 2n7r and (2n + l)7r, 
respectively. The blue dashed line represents the behaviour 
at K = and is plotted for reference. The other parameters 
are Vl = O.O57, A = 0, i/jl = 0. 

To conclude this section, one of the main features as- 
sociated with photon-mediated atom-atom interaction is 
that the intensity-intensity correlation exhibits an en- 
hanced or suppressed probability to measure a second 
photon as a function of the time after the first detection. 
This behaviour is due to interference between the various 
paths of multiple scattering, and can be interpreted as a 
combined photonic-atomic excitation which is stored in- 
side the system. In view of the interpretation by [13, HI] , 
one can say that for a transient time the system develops 
and stores entanglement and correlations, determined by 
the strength of the interaction k, until the atoms finally 
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dissipate the excitation into free space. In the future, 
it would be interesting to consider these dynamics in a 
quantum jump picture [s^l- This could open the possi- 
bility of implementing schemes for entangling atoms in 
this kind of set-up, as proposed in [s^]. 

VI. CONCLUSIONS 

We have studied the photonic properties of two atoms 
which are coupled by radiation via an optical element 
such as, e.g., a lens or an optical fiber focussing a rele- 
vant fraction k of electromagnetic field modes from one 
atom to the other. Signatures of multiple scattering of 
photons between the atoms are observed in the first-order 
and second-order coherence of the scattered light. These 
features show that the presence of the second atom sub- 
stantially modifies the radiative properties of the first 
one, even when the atoms are separated by a distance d 
much larger than the light wavelength A. 

The efficiency of the interaction, which for two atoms 
in free space scales with X/d, is determined by n when 
an optical system mediates the coupling. The atom-atom 
distance d plays a new role, separating two regimes where 
the delay of the interaction t — d/cis smaller or larger 
than the atomic decay time I/7. In these two regimes 
7T < 1 and 7T > 1, the coupled two-atom system shows 
characteristics of a single- or multi-mode resonator, re- 
spectively, with mirrors of low reflectivity k, and band- 
width 7. 

In this article we considered the limit in which the 
atoms are weakly driven by the laser, and we neglected 
the effect of atomic motion. It should be remarked that 
localization of the atoms within a wavelength of the 
scattered light is a relevant requirement for observing 
the interference effects arising from multiple scattering. 
In fact, atomic motion gives rise to a dephasing in the 
signal, which can be interpreted as which-way informa- 
tion imprinted by the photon recoil on the scattering 
atom |32|, |39|, H^]. Nevertheless, when the recoil of the 
atom in each photon scattering event is taken into ac- 
count, correlations between the atomic motion and the 
light are established [4l[. In particular, mechanical ef- 
fects between the distant atoms arise, which are medi- 
ated and retarded by the optical coupling [l^ |42| • It is 
interesting to consider whether such effects may lead to 
novel collective behaviour of atomic center-of-mass and 
photonic variables, in analogy to collective dynamics pre- 
dicted for cold atoms inside resonators I43l. 
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APPENDIX A 

We formally integrate Eqs. (j9bp and ((9c)) using the 
initial condition h'f\Q) = 0. Inserting the result into 
Eq. ^ yields 

t 

h^pit) = - ^ gj j di'e'("«-"^)(*-*')&ii)(t') (Al) 
e.=ij.,p Q 

t 

- H 5p / dt'e'("°-"'')(*-*')e*''-("-i-'-^)6(2) {t') , 

" 

where the equation for be (t) is found by swapping the 
indexes 1 and 2 in Eq. (jAip . We see that the differential 
equation for the probability amplitude bi^\t) depends 
linearly on the probability amplitude be {t') for the 
excited state of atom 2. Such dependence is due to the 
common modes p which mediate the interaction between 
the two atoms. In absence of the second atom, the 
equation reduces to the well known equations describing 
the radiative decay of a two- level atom in free space [27[ . 

The first term on the right-hand side (RHS) of Eq. (|XT|) 
can be rewritten in a compact way by converting the 
sum over the modes into an integral. Using the Wigner- 
Weisskopf approximation one obtains [3| 

h^e\t)^-lbP{t) (A2) 
t 

- j dt'bP{t')J29y"''^"''"^^e'^'""'^'^^'~*'K 
P 

where 7 = D^lu^/ (STreohc^) is the free-space decay 
rate . The second term on the RHS of Eq. corre- 
sponds to the sum of the modes mapped from one atom 
into the other by the optical setup. Let us consider a lens 
between the atoms collecting a solid angle S^Iq of modes 
with aperture ^o- Converting the sum over the modes 
into an integral, the sum over the modes p in Eq. (|A2p 
can be rewritten as 

(^2gikp-(ri-r2)gi(wo-'^p)(t-t') _^ 

P 

00 

2(27r)3eofc3 J 


X j dr!e''' ("-^--=) [d^ - • (A3) 
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Since optics compensate for the phase difference between 
the various modes, we take 



j,ik(ri-r2) 



(A4) 



where r is defined in Eq. ([T]), and is the time a pho- 
ton emitted inside the sohd angle S^Iq needs to cover 
the distance between one atom and the other via opti- 
cal setup. Using Eq. (jA3|l into Eq. (jA2|) . we can now 
make the Wigner-Weisskopf approximation and obtain 
Eq. dlOD. 



As we are interested in evaluating the scattering pro- 
cesses for long times, t —^ oo, we neglect the term kr in 
the upper bound of the integral and take the Heaviside 
fimction to be one. The first term inside the parentheses 
on the RHS of Eq.(|B6| gives 



dxe 



-iA„ 



g-iA,t/22^5(t)(^^) 



(B8) 



In order to evaluate the term (jB7|l we use the relations 
Eq. jHll and 



APPENDIX B 

We consider the terms of Eq. (I??)) . Let us introduce 
the simple relation 

t 

1= f dt'S^^~^-'>''Hk{t',u;^) 



t-kr 



k\ 



''^^''^e{t^ kr) / dxe'('^^-"^)^x'=Gfe(aa;), 



(Bl) 



with a ^ ^ + i6fi and where we used Eq. ([^ . Usin g th e 
definition of the confluent hypergeometric function [47[, 



^ , , , a a(a + 1) 



we rewrite Eq. pT|) as 



or equivalently 



Gk{s) 



k ^ — s 



e 



(B2) 



1 - e-^ if fc = 0; 



(B4) 



; x-'e 



j! i!e-"' ^ {aty 







1=0 



in Eq. (|B7p . which then reads 

k 



Skit) = 



a-' 



\ 1=0 



(« + iAp)'t' 

n 



In particular, 



02k 



k=0 

where 



K'^''+'S2k+i " 2 



J2K"'Sk±{-l)''K"'Sk 



Lk=0 



J2{±^)''K"'Sk 



k=0 



(B9) 
(BIO) 



with K' = — K^e'"^'^/(^ -I- ii5^), and where we used the 
Cauchy-Product for absolute convergent series, thereby 
neglecting the vanishing exponentials as we consider the 
long time limit. 

In the long-time limit, using that limj^+oo = S{x) 

we finally arrive to the relation in Eq. ()3ip . 



whereby 



APPENDIX C 



n— ^ ^ 

We use Eq. dH]) in Eq. (|Bl]), and obtain 



■ t-kr 



' Q{t - fcr) J dxe'("--'^")^ - Sk (96) 



with 



^^^-W^^r / (B7) 



Starting from Eq. (H^l) we calculate the second or- 
der correlation function of two atoms, which are weakly 
driven by the laser and both scatter towards the detector. 
In the reference frame rotating at the laser frequency, as- 
suming that the initial state is the atomic ground state, 
we rewrite Eq. ([44]) as 



e(2)(xi,<;X2,< + t') = (CI) 

II (ai + a2e}nU{t'){cT^ + a^e^'^'Wm, g, 0) |p , 

where the correlation function is evaluated at lowest or- 
der in perturbation theory in the atom-photon interac- 
tions. The operator U{t) is the total evolution operator, 

U{t) = exp{-iH't/h) , 
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with H' given in Eq. ([5]), which is to be expanded in Inserting the exphcit value of the coefRcients, Eq. (I25p . 
power series of the interactions V^mf and V^. At lowest, after partially integrating B{t) reads 
non-vanishing order, Eq. (jCip can be rewritten as 



= \{A^{t) + A2{t)e''^') {A^{t') + A^ity'"^') + B{t) 



(C7) 



(C2) 



|'l-,2ikL ri ^ ^2ikL r2 



where the coefficients Aj and B are the transition ampli- 
tudes 

Ai(t) = (e,g,0|C/(0l5,5,0) = c«(i), (C3) where 
A^it) = (5,e,0|C/(t)|.9,.g,0) -c(2)(t), (C4) 
B(t) = (e,e,0|C/(t)|.9,.g,0), (C5) 

where ci\t) is given in Eq. ((25|) . while the probability i^j, — 
amplitudes Cji are defined as 





= (e,g,0|C/(t)|e,g,0) 








= e-'^%\}\t) with 




= 1, 


Cl2 


= (.g,e,0|C/(t)|e,g,0) 












= 1, 


C21 


= (e,g,0|{7(t)|.9,e,0) 








= e-~'^%\}\t) with 


&(^)(0) 


= 1, 


C22 


= (.g,e,0|C/(t)|.9,e,0) 










6(^)(0) 


= 1, 



with bi^\t) given in Eq. p^ . We evaluate the coefficient 
in second-order perturbation theory, hence obtain- 
ing 



k 



7 — KC 

2 ^ 



^,Q(t - fcr)(r(fc + 1, a{t - fcr)) 



+ (-l)'=e-2"(*-Mr(fc + 1, -Q(f - fcT))| 



with a = + iA and r(fc, a) the generalized gamma 
function [43. In the long-time limit the second term in 
the curly brackets is negligible, and we find 



B{t) 



2(^+iA)2 l-if2 



(C8) 



2gikL (ri+r2) _ ^^2ikL ri _^ ^2ikL r2-J 



B{t) = -in / dt'e-(2'^+^)(*-*') 



while the form Aj (t) valid in the long-time limit is given 
^ by Eq. ([30| . Inserting the explicit value of the coefficients 

U^'^"^ci^\t') + e"'^ ''ic(2)(t')j . (C6) in Eq. ^ we finally obtain Eq. gb^. 
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